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Abstract 

The rapidity dependence of the single- and double- neutron to proton ratios of nucleon emission 
from isospin-asymmetric but mass-symmetric reactions Zr+Ru and Ru+Zr at energy range 100 ~ 
800 A MeV and impact parameter range ~ 8 fm is investigated. The reaction system with isospin- 
asymmetry and mass-symmetry has the advantage of simultaneously showing up the dependence on 
the symmetry energy and the degree of the isospin equilibrium. We find that the beam energy- and 
the impact parameter dependence of the slope parameter of the double neutron to proton ratio (Fr>) 
as function of rapidity are quite sensitive to the density dependence of symmetry energy, especially 
at energies Ef, ~ 400 A MeV and reduced impact parameters around 0.5. Here the symmetry 
energy effect on the Fd is enhanced, as compared to the single neutron to proton ratio. The degree 
of the equilibrium with respect to isospin (isospin mixing) in terms of the Fr> is addressed and its 
dependence on the symmetry energy is also discussed. 

PACS numbers: 24.10.Lx, 25.75.Dw, 25.75.-q 
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High energy collision of heavy ions provides a unique tool to study the unbroken prop- 
erties of dense nuclear matter at finite temperature, outside of the reach of conventioned 
nuclear physics . The FOPI Collaboration has recently performed the mixing experi- 
ment by using four mass 96+96 systems 44RU+44RU, ^Zr+^Zr, ^Ru+^Zr, and ^Zr+^Ru 
to investigate the degrees of the isospin-equilibrium in the intermediate energy heavy ion 



collisions (HICs) 



Aration^is not reached by 
). The advantage 




I. It was found that complete isospin equi 
eo« the emitted proton nUmb e r counted to r apidi ty J 
of taking mixed reactions is that the isospin asymmetry shows up more obviously because 
the effects of the iso-scalar part of the equation of state (EoS) largely cancel. This advan- 
tage could also be used for probing the density dependence of the symmetry energy. So far 
the density dependence of the nuclear symmetry energy, E sym (u), with u = p/p being the 
reduced nuclear density, is still largely uncertain, although a lot of theoretical and experi- 



mental efforts have been undertaken in the past decade |8, 



Recently, the symmetry 



energy has been extracted to E, 



sym 



U 



31.6m 105 at subnormal densities by comparing with 



experimental data |TJ| when a free nucleon-nucleon (NN) elastic cross section has been used 
|l2j ]. while E Lymju ) — 31.6m ' 69 after considering partly medium modifications of the nucleon 
transport |13l Il4j ] were subtracted . The symmetry energy can be investigated by the double 
neutron to proton and 7r~ /tt + ratios. The neutron-rich 124,132 Sn + 124 Sn system and more 
isosp in-s ymm etric 112 Sn + 112 Sn system serve as robust probes of the nuclear symmetry en- 



ergy 



HE 



A very attractive feature of taking double ratios is the fact that the systematic 
errors can be strongly reduced, while the sensitivity of the symmetry energy is close to the 
single neutron to proton and tt~/tt + ratios. 

In the present work, we take two sets of isospin-asymmetric but mass-symmetric reactions, 
Zr+Ru and Ru+Zr, to study the symmetry energy effect on the rapidity dependence of both, 
single as well as double unbound neutron to proton ratios. We take advantage of the mixing 
reactions to study the influence of the isospin dependent part of the EoS on the equilibration 
process. We show that the slope of the double neutron to proton ratio as distributed in the 
rapidity space is — for mixing reactions — closely related to degree of isospin equilibration. 



re UrQMD transport model has been updated for simulating intermediate energy HICs 




18, 



19 



y 



2l|. It is also used for the calculations presented in this work. A soft 
equation of state (S-EoS) is adopted 2qJ (except otherwise stated). We consider here two 



sets of symmetry energy parameterizations: E, 



sym 



u 



34m 7 with 7 = 0.5 (soft) and 1.5 
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(hard). A non-relativistic medium modification of the nucleon-nucleon elastic cross section 
is also considered here (please see Ref. j2jj for details). The method to construct the 



freeze-out clusters is the same as presented in our previous work ljj, [21 1 . 

We calculate the rapidity {Ycm = VcmlVbeam hi the center-of-mass system) dependence of 
unbound neutrons and protons in Zr+Ru and Ru+Zr reactions at the beam energies E b = 
100, 200, 400, 600, and 800 A MeV and at the reduced impact parameters b/b = 0, 0.25, 
0.5, and 0.75, where bo = R pro j + Rtarg is the maximum impact parameter leading to nuclear 
reactions by adopting 'sharp-cutoff' approximation. The calculated rapidity dependence of 
the neutron to proton "n/p" ratios for Zr+Ru and Ru+Zr reactions are shown in Fig.ffl(a)- 
(d) for E b = 100 A MeV, b/b = 0.5 in plot (a); E b = 400 A MeV, b/b = 0.5 in (b); E b = 400 
A MeV, b/b = in (c); and E b = 800 A MeV, b/b = in (d). If equilibrium with respect 
to the isospin degree of freedom was reached, the rapidity distributions of the neutron to 
proton ratios for the two different mixing reactions should coincide with each other. Figs. 
E](a) and (b) show clearly that for b/bo = 0.5 equilibrium (with respect to isospin degree of 
freedom) is not reached, because the centroid of the rapidity dependence of the neutron to 
proton ratio is not located at mid-rapidity |5]. For Zr+Ru reactions, the neutron to proton 
ratio distribution peaks at the positive rapidity, while for the inverse Ru+Zr reactions, the 
distribution peaks at the negative rapidity. The rapidity dependence of the neutron to proton 
ratio is symmetric with respect to Ycm = for central reactions at high beam energy (Fig. 
12(d)), where the equilibration is nearly reached. It is interesting to see that the statistics 
in the UrQMD calculations is perfectly guaranteed since the calculated results for Zr+Ru 
and Ru+Zr reactions are symmetrically distributed very well with respect to Ycm = 0, as 
one expects. This is also a fundamental but important checkup in the experimental side. As 
seen for reactions at E b = 400 A MeV and b/bo = in Fig. ^ (c), the rapidity distribution 
of the neutron to proton ratio almost coincides for reactions Zr+Ru and Ru+Zr when the 
soft-symmetry energy is chosen, but it deviates strongly when the stiff-symmetry energy is 
chosen. This shows the closely relationship between the effect of the symmetry energy and 
the equilibration in HICs at intermediate energies. 

Figs. ^ (e) and (f) present the corresponding rapidity distributions of double neutron 
to proton ratio "(n/p)z r Ru/( n /p)RuZr" with the same conditions as Figs. [T] (a) and (d), 
respectively. From Fig. ^ (e) it is seen that the double ratio increases almost linearly from 
target {Ycm = —1) to projectile {Ycm = 1) rapidity region for the case of E b =100 A MeV and 
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b/bo = 0.5. The two curves corresponding to soft and stiff symmetry energies cross at the 
point of rapidity equal to zero, with positive but different slopes. And obviously, the slope 
for 7 = 0.5 case is larger than that for 7 = 1.5 case. It is clear that for isospin asymmetric 
and mass symmetric reactions Zr+Ru and Ru+Zr, if the equilibrium with respect to isospin 
degree of freedom is totally reached, the double ratio should be unity at all rapidities. 
The deviation of (n/p) zr+Ru/ \ n l p) Ru+Zr from unit thus implies the deviation from isospin 
equilibration. Fig. ^ (e) shows clearly that the isospin equilibration is not reached and the 
choice of symmetry energy influences the isospin equilibrium process for this case. While for 
the reactions at Eb = 800 A MeV and b/bo = presented in Fig. [T] (f), the slope is ~ and 
becomes horizontal in large rapidity region, which means that the isospin equilibration can 
be roughly reached. Thus the double neutron to proton ratio as function of rapidity taking 
from the mixing reaction system can provide both the information of density dependence of 
the symmetry energy and the degree of isospin equilibrium reached in the reaction. 

Further, we study the beam energy and the centrality dependence of the slope parameters 
of the double neutron to proton ratio over the rapidity, namely the Fd, which is obtained by 
linearly fitting the rapidity distribution of the double ratio over \y£}\ < 1. Fig. (a) shows 
the excitation function of Fd for impact parameters b/b = and 0.5 and for symmetry 
energies 7 = 0.5 and 1.5, respectively. The reactions with larger impact parameters show 
larger Fd values, which implies that less equilibration is reached. The Fd values first rise and 
then drop slightly with increasing beam energy. It shows a broad maximum at Eb ~ 200—400 
A MeV for intermediate impact parameters, while central collisions peak at 200 A MeV. The 
influence of the symmetry energy on the Fd is stronger for b/bo — 0.5 than for b/bo — 0. The 
isoscalar part of the EoS affects the equilibrium process of the colliding system as well j^. 
Fig. |21 (a) shows Fd at Eb = 400 A MeV and b/bo = 0.5 using a soft EoS with momentum 
dependence (SM-EoS, open dots). It is seen that the difference between the S-EoS and the 
SM-EoS is very small. This indicates that the effect of the isoscalar part of the EoS on Fd 
is largely cancelled, while the effect of the isovector part on Fd is nicely preserved. 

Fig. El (a) also shows that — for central collisions — Fd as calculated with a stiff symmetry 
energy is larger than with a soft one. For b/bo — 0.5, the relative magnitude of Fd as 
calculated with the stiff symmetry energy becomes smaller than that with the soft symmetry 
energy. The centrality dependence (with b/b Q = 0, 0.25, 0.5, and 0.75) of F D at beam energy 
400 A MeV is presented in Fig. E](b) . This shows that both curves of Fd, as calculated with 
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FIG. 1: (Color online) (a)-(d): Rapidity dependence of neutron to proton (n/p) ratios with sym- 
metry energies of 7 = 0.5 and 1.5. The Zr+Ru and Ru+Zr reactions with E b = 100 A MeV, 
b/b = 0.5 in plot (a), E b = 400 A MeV, b/b = 0.5 in (b), E b = 400 A MeV, b/b = in (c), and 
E b = 800 A MeV, b/bo = in (d) are chosen, (e)-(f): the corresponding double neutron to proton 
ratios ((n/p)z r Ru/ ' {n / P) RuZr) for (a) and (d) cases (see text). 
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soft and stiff symmetry energies, increase monotoneously with impact parameter. The Fd 
value obtained with the stiff symmetry energy is larger than that with the soft symmetry 
energy, when b/bo < 0.2. And they reverse when b/b > 0.2. The reversal of the relative 
magnitudes of Fd with stiff or soft symmetry energy at small and large impact parameters 
is due to the gradual change of the mechanism of particle emission U. Central collisions 
are more violent, hence the compression is stronger than for the collisions with large impact 
parameters. Correspondingly, the degree of isospin equilibration is both higher and reached 
more quickly. Nucleons emitted from central collisions stem mainly from the mid-rapidity 
region, where matter is much more compressed than at the projectile/ tar get region. During 
the compression phase, the stiff symmetry potential provides stronger repulsion (attraction) 
than the soft one for neutrons (protons) in the neutron-rich matter. Therefore, the center 
of collisions is more strongly compressed. However, higher pressure caused by the stiff 
symmetry potential leads to faster emission of particles. Hence, fewer two-body collisions 
occur and the equilibration decreases. For peripheral reaction, the average compression is 
relatively weak, Nucleon emission takes place mainly at projectile/ tar get rapidity, where 
the nuclear density is low. The low-density nuclear surface thus plays an important role for 
nucleon emission. At low densities, the stiff symmetry potential provides weaker repulsion 
(attraction) on neutrons (protons) than the soft one in the neutron-rich nuclear medium. 
Higher densities result in the projectile/target region, leading to an enhanced number of 
two-body collisions. The smaller Fd value is evident with the stiff symmetry energy. 

The effect of the symmetry energy on Fd becomes the most pronounced at Ef, ~ 400 A 
MeV and at large impact parameters (Fig. EJ). Compare the isospin effects in Fig. [2(b) and 
Fig. |2]for reactions at =400 A MeV and b/b = 0.5: The double neutron to proton ratio 
exhibits an enhanced sensitivity (~ 30%) to the density dependence of the symmetry energy 
as compared to the single neutron to proton ratio (~ 10%). Thus, the double neutron to 
proton ratio (taken from mass symmetric, but isospin asymmetric mixing systems at energies 
around 400 A MeV and at semi-peripheral impact parameters) is a prominent probe of the 
density dependence of the symmetry energy. 

In summary, we have investigated the rapidity dependence of the single- and double 
neutron to proton ratios as taken from isospin-asymmetric but mass-symmetric systems 
Zr+Ru and Ru+Zr. The beam energy- and the centrality dependence of the slope parameters 
of the double neutron to proton ratio (Fd) as function of rapidity are presented. We find 
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FIG. 2: (Color online) (a): Beam energy dependence of the slope parameter Fjj with symmetry 
energies with 7 = 0.5 and 1.5, for impact parameters b/bo = and 0.5. The Fd values at Ef, = 400 
A MeV and at b/bo = 0.5 with SM-EoS are shown with open dots, (see text) (b): The centrality 
dependence of F D at E b = 400 A MeV. 

that both, the rapidity distribution of the single neutron to proton ratio and the Fd are 
sensitive to the density dependence of symmetry energy. The study of the beam energy- 
and the centrality dependence of the Fd values shows that the (semi-) peripheral HICs at 
energies around ~ 400 A MeV are most suited for probing the density dependence of 



S 



the symmetry energy. Here the symmetry energy effect is enhanced most strongly as is 
observable in the rapidity dependence of the double neutron to proton ratio, as compared 
to the single ratio. 
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